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Abstract 

A general framework for studying a large class of cosmological solutions of the low-energy 
limit of type II string theory and of M-theory, with non-trivial Ramond form fields excited, is 
presented. The framework is applicable to spacetimes decomposable into a set of flat or, more 
generally, maximally symmetric spatial subspaces, with multiple non-trivial form fields spanning 
one or more of the subspaces. It is shown that the corresponding low-energy equations of motion 
are equivalent to those describing a particle moving in a moduli space consisting of the scale 
factors of the subspaces together with the dilaton. The choice of which form fields are excited 
controls the potential term in the particle equations. Two classes of exact solutions are given, 
those corresponding to exciting only a single form and those with multiple forms excited which 
correspond to Toda theories. Although typically these solutions begin or end in a curvature 
singularity, there is a subclass with positive spatial curvature which appears to be singularity 
free. Elements of this class arc directly related to certain black p-brane solutions. 
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1 Introduction 



An important constraint on string theory or any generalization of string theory, such as M-theory, 
is that it should be compatible with the standard model of early universe cosmology. The standard 
approach to string cosmology has been to consider an epoch when the evolution of the universe 
is described by the string low-energy effective action. One encounters a number of the familiar 
problems in matching this description to standard cosmology pj|. In particular, it is necessary 
to invoke some mechanism for stabilizing the massless string moduli such as the dilaton. While 
understanding these issues may in fact be central to understanding string cosmology Q , the usual 
approach has been at least to start by studying the simple cosmological solutions of the pure 
effective supergravity theory. 

In the past, most focus has been on heterotic strings as the best model of low-energy particle 
physics. Consequently, study of cosmological solutions has concentrated on the dynamics of the 
dilaton and the compactified space in cosmological solutions of the effective theory. However, with 
the discovery of string dualities || and the existence of D-brane states |J, the nature of string 
theory has changed dramatically. Strong-weak coupling duality symmetries connect each of the 
five consistent super symmetric string theories together with eleven-dimensional supergravity ||. 
As a consequence, type II and eleven-dimensional supergravities may now be directly relevant to 
low-energy particle physics and cosmology ||. Both theories contain form fields, namely a three 
form in eleven-dimensions and Ramond-Ramond (RR) forms of various degrees in type II theories. 
The importance of considering the excitation of these fields has been stressed by the discovery of 
objects carrying RR form-field charge, the D-branes of open string theory, and the central role 
these states have played in understanding the spacetime structure of black holes |7j. Given this 
change of perspective, it clearly becomes important to study the cosmological solutions of type II 
and eleven dimensional supergravity with non-trivial form fields excited. In a recent letter ||, we 
gave the first study of a class of such solutions, giving specific examples with one or more forms 
excited. The purpose of this paper is to provide a more detailed and complete description of this 
class of solutions. 

The gravitational part of our Ansatz is based on the standard "rolling radii" solutions of Mueller 
which are a generalization the classical Kasner solutions of Kaluza-Klein gravity fllofl . In the 
original Mueller solutions spacetime is divided into a set of maximally symmetric spatial subspaces, 
each with its own time-dependent scale factor. (One subspace should be three-dimensional to 
describe our observed spacetime.) The rate of expansion or contraction of the subspaces is controlled 
by a power law and, typically, the universe either starts or ends in a curvature singularity. Including 
a dilaton in ten dimensions does not significantly alter the solution. This is not surprising given that 
one may consider the dilaton as the radius of a compact eleventh dimension. As was first pointed 



out in [11 1, in general such an evolution describes a motion on a moduli space with coordinates 
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given by the values of the scale factors together with the dilaton. (We should point out that 
this is not quite the conventional notion of a string moduli space. For instance, it includes a scale 
factor "modulus" for our observed four-dimensional spacetime, and does not include all the possible 
moduli of, for example, toroidal compactification, since it assumes the torus is diagonal.) 

We will adopt this rolling radii Ansatz for the metric and the dilaton. However, the point of 
this paper is to study solutions of type II theories and M-theory by turning on non-trivial Ramond- 
Ramond form-fields. Can we understand how these will modify the moduli space description? Such 
solutions were first considered in the context of eleven-dimensional supergravity by Freund and 
Rubin [12|. Solutions with a non-trivial Neveu-Schwarz two-form field, usually in four dimensions, 
have also been discussed by various authors ||, 0, [jj. Following Freund ||, if we are 

to preserve the structure of a spacetime split into a set of maximally symmetric subspaces, the 
form fields must be chosen to span completely one or more of the subspaces. They may or may 
not in addition have a component in the time direction. Our central result will be to show that 
the effect of such a configuration is to introduce new exponential potential terms into the moduli 
space equations. In general the new equations describing the evolution of the dilaton and the scale- 
factors can only be solved exactly in certain special cases. We presented some examples of these 
special cases in a previous letter ||. (One example appeared at almost the same time in paper by 
Kaloper A subsequent paper by Lii et al p9[ , gave a broad class of further exact solutions. 

Recently, cosmological solutions with Ramond forms obtained from black hole solutions have also 



been studied [16, 20 1. The present paper provides a general framework for analyzing cosmologies 
with Ramond forms, the Neveu-Schwarz form and the cosmological constant of type II theories as 
well as those in eleven dimensions. 

For all the cases we can solve exactly, our general result is that the moduli-space potential which 
arises when non-trivial form-fields are excited, significantly effects the structure of the solution. 
When all the subspaces are flat, the potential is operative over a particular, finite part of the 
evolution; at the extremes, the solution becomes pure geodesic motion and we return to the simple 
Kaluza-Klein Kasner-type solutions with some subspaces expanding, some contracting. (There is 
consequently always either an initial or a final curvature singularity.) Thus the effect of the form 
fields is to interpolate between two different Kaluza-Klein solutions. As we have mentioned, the 
expansion or contraction in the Kaluza-Klein solutions is controlled by a power law. In general, 
for solutions with an initial singularity, one finds that the rate of expansion is always sub-luminal 
and so there is no inflation. On the other hand, those solutions with a final singularity, just as in 



the pre-big-bang models |21], may exhibit super inflation but are unphysical as stands because the 
inflation ends in a curvature singularity. 

The solution can be very different, however, when we allow for curved subspaces. The effect 
of the curvature is to introduce new terms into the moduli space potential, with the particular 
new feature that, in the case of a spherical subspace, the new term is always negative. This can 
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significantly change the singularity structure, in some cases giving solutions where the curvature 
always remains finite. This suggest the interesting possibility that there may exist inflating solutions 
which are not forced to end in a singularity. 

It is worth noting that there is actually a close connection between cosmological solutions and 
those describing p-branes, as was pointed out for the first time in ref. || and recently analyzed, 
starting from black hole solutions, in [20|. In the single p-brane solution J5|, all the fields depend 
on the radial direction transverse to the brane. If the brane is elementary the form field has a 
component in the radial direction. If it is solitonic it spans a sphere around the brane, with no 
radial component. Since our solutions depend on time alone, time must play the role of the radial 
coordinate. Thus these solutions should correspond to purely spacelike p-branes, spatial (p + 1)- 
dimensional objects which appear at one instant in time. Similarly, the two possible orientations of 
the form field in our solution, with a time component or without, correspond to the elementary and 
solitonic branes respectively. There is really one subtlety in this interpretation, at least when all 
the subspaces are flat. A standard p-brane solution in D dimensions has a {D — p — l)-dimensional 
transverse space, composed of a radial direction and a {D — p — 2)-sphere, whereas, in the cosmo- 
logical solutions, time is the only transverse direction. However it is well known that the number of 



effective transverse directions can be reduced by stacking arrays of parallel p-branes together [22 



In the limit of smooth distribution of branes the fields no longer depend on the transverse coordinate 
along the direction of the array but only on the coordinates transverse to both the brane and the 
array. Thus we see that our cosmological solutions should really correspond to a distribution of par- 
allel, spacelike branes evenly spread throughout space, though again appearing at only one instant 
in time. Solutions with multiple form fields excited then correspond to distributions of intersecting 



spacelike brane [23]. It is known that p-branes can be considered as interpolating solutions between 
different purely gravitational backgrounds. Therefore, the result that our cosmological solutions 
interpolate between pure Kaluza-Klein solutions supports the above analogy. The analogy breaks 
down however with regard to super symmetry. Unlike the case of conventional p-brane solutions, we 
found no cosmological solutions which preserved any supersymmetries. We will see that the actual 
form of the solutions is closer to the black p-brane |p4], || and other non-extremal brane solutions 

If one allows non-flat subspaces the p-brane analogy can become even clearer. Behrndt and 
Forste [|l(| and, more recently, Poppe and Schwager [20| have obtained cosmological solutions 
directly from black p-brane solutions by going to a regime where roles of the time and radial 
coordinate are exchanged. Such solutions correspond in our formalism to allowing a single spherical 
subspace. This subspace together with the time direction then constitute the space transverse to a 
single p-brane. This raises the possibility of finding new cosmological solutions of this form which 
correspond directly to single black p-brane states. 

The organization of the paper is as follows. In section two we give the general structure of the 
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solutions. Starting with the form of the effective action, we give the metric Ansatz and discuss the 
possible elementary- and solitonic-type Ansatze for the form fields. We then describe the action in 
moduli space which reproduces the spacetime equations of motion, and discuss its form. In section 
three, we give the exact general solution corresponding to exciting a single form field with flat 
subspaces, and discuss its structure as an interpolating solution between two pure Kaluza-Klein 
solutions, together with an explicit example. In section four, we give a class of exact solutions with 
multiple forms excited. These are solvable because they correspond to Toda models. Again, giving 
explicit examples in two distinct cases with flat subspaces, we discuss the interpolating form of 
these solutions. In section five, we explain how to incorporate the case of curved spatial subspaces 
into our formalism. A class of solutions appears to be related to black p-branes. We give, as an 



example, the solution of Behrndt and Forste [16], and discuss the possibility that those solutions 



which are related to black p-branes do not begin or end in a curvature singularity. We present some 
very brief concluding remarks in section six. 

2 General framework 

In this section we present a general framework for finding cosmological solutions with non-trivial 
form fields. While here the discussion will be as general as possible, we briefly discuss the same anal- 
ysis in terms of simple, illustrative example in section 3.2. The starting point for our investigation 
is the following effective action 



(i) 



2{8 r + 1)! 

with the L>-dimensional string frame metric <jmn, the dilaton </>, the NS 2-form H and a number of 
RR 5 r — forms F r = dA r . We also allow for a cosmological constant A which appears in the massive 
extension of IIA supergravity [26|. Assuming this origin, it is restricted to be positive, A > 0. Note 



also that this cosmological constant behaves much like a RR form in that is does not couple to 
the dilaton in the string frame. In this respect it is quite different from a "normal" cosmological 
constant appearing in an effective action of noncritical dimension which, like the NS field, couples 
to exp(— 2(f>). 

The above action can account for a wide range of cosmological solutions in type II theories 
(where we usually have D = 10 in mind) and, if the dilaton is set to zero, in D = 11 supergravity. 
For simplicity we have kept only the kinetic terms for the form fields. In general, in both type II 
theories and eleven dimensional supergravity there are additional terms involving the coupling of 
form fields. We shall assume throughout this paper that for the configurations with which we are 
concerned, these terms do not contribute to the equations of motion. This will always be true if 
we only excite one orientation (in the sense given below) of a single form. It is not necessarily true 
when multiple forms are excited. 



4 



In order to give a physical description of our solutions, we prefer to work in the canonical 
Einstein frame metric qmn which is related to the string frame metric by a conformal rescaling 
9mn = exp(— 4/(_D — 2)4>)gMN- The corresponding action reads 



S = I d D x 



R - (d6) 2 - V - e ~ a ( s r)<t> F 2 _ A e ~ a A^ 

D-2 {9) ^2(<J r + l)! 



(2) 



where the NS field H has been included in the sum over F r . It is distinguished from the other 
forms by the dilaton couplings a(S r ) given by 

NS 2 - form 
RR <L— form 



= < ^-2^-2) : : • ( 3 ) 



£1-2 

The coupling a\ for the cosmological constant 

a A = (4) 

A D-2 v ' 

equals the negative of the one for a RR (D — l)-form. This reflects the fact that it couples to the 
8-brane solution [^J of IIA supergravity. 

The type of solutions we consider are characterized by a space split into n flat subspaces, each of 
them characterized by a scale factor on. We concentrate on flat subspaces since this covers already 
a wide range of solutions and seems to be appropriate for our main purpose, namely to investigate 
the effect of the forms. The effect of curved maximally symmetric subspaces can, however, be easily 
incorporated into the formalism to be developed in this section. We will discuss this in section six. 
In the flat case, the corresponding metric is given by 

n-i 

ds 2 = -N 2 (r)dT 2 + J2 e 2ai(r) dx- (5) 



i=0 



where dxf is the measure of a dj-dimensional flat maximally symmetric subspace and J27=o ^» = 
D — 1. For solutions with this structure, the dilaton should depend on time only, 4> = <M r )- This 
Kaluza-Klein-type Ansatz is about the simplest allowing for the cosmologically key properties of 
homogeneity and isotropy as well as for an "external" and "internal" space. 

For a ten dimensional theory, the simplest choice is to split up the space into two subspaces 
(n = 2) with do = 3 and d\ = 6. Then the do = 3 part could be interpreted as the spatial 
part of "our" 4-dimensional space-time with an evolution described by the scale factor olq. The 
other six dimensions would form an internal space with a modulus a\. Clearly, one is free to 
split these six dimensions even further and to consider, for example, a situation with n = 3 and 
(do, d%, ^2) = (3, 4, 2). From the cosmological point of view one might even allow for a further split 
of the 3-dimensional space such as (<io>^i>^2) = (1)2,6) as long as the 1 + 2-dimensional space 
approaches a 3-dimensional isotropic space at late times. 
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As discussed in the introduction, this Kaluza-Klein-type Ansatz is closely related to the form 
of p-brane solutions in supergravity theories. Consider the simple case of two subspaces with 
<^o = P + 1 an d d\ = D — p — 2. For the standard p-brane solution, spacetime is split into a radial 
direction, transverse to the brane, upon which all the fields depend, a flat subspace describing 
the (p + l)-dimensional worldvolume of the brane and hence including a timelike direction, and a 
spherical subspace surrounding the brane. In our cosmological solution the fields depend only on the 
timelike direction and so time is the analog of the radial coordinate in the p-brane solutions. The do 
subspace then becomes the analog of the brane worldvolume, though it is now completely spacelike, 
while the d\ subspace becomes the sphere surrounding the brane. As described in the introduction, 
since we usually take the d\ space to be flat rather than spherical, probably a more exact analogy 
is with a distribution of spacelike p-branes, rather than a single brane. More complicated cases, 
where the space is split into a number of subspaces, are analogous to multiple intersecting p-brane 
solutions [p3[| . 

The choice of subspaces is important in fixing the possible structure of the antisymmetric tensors 
fields. The symmetry of the above metric allows for two different Ansatze for the form fields which 
we call "elementary" and "solitonic" in analogy to the two types of p-brane solutions. They are 
characterized by the following nonvanishing components of the field strength. 

• elementary : if J2i d% = S r for some of the spatial subspaces i we may set 

(i ? r)o Att ... Wr = A(a)fUr)e^ Sr , A r (a) = e" 2 ^** (6) 

where fMi...fj,g r refer to the coordinates of these subspaces, f r (T~) is an arbitrary function to be 
fixed by the form field equation of motion, and the prime denotes the derivative with respect 
to t. With raised indices, the symbol e w '" Wr takes the values or 1 and is completely 
antisymmetric on all 5 r indices. Note that the sum over i in the exponent runs only over 
those subspaces which are spanned by the form. 

An elementary form can therefore extend over one or more of the subspaces only if its degree 
matches the total dimension of these spaces. Consider for example the RR three-form of type IIA 
supergravity. If the space is split up as (do,di) = (3,6) it fits into the 3-dimensional subspace and 
the above general Ansatz specializes to Fq^^^ = exp(— 6ao)/'(T)e AllA12M3 where ^i,//2,A*3 refer to 
the coordinates of this subspace. Taking (do,di,d2) = (1,2,6) allows for a similar Ansatz but the 
form now extends over the 1 + 2-dimensional space and the exponential in the Ansatz changes to 
exp(— 2ao — 4ai). 

• solitonic : if J2i d% = 5 r + 1 for some of the spatial subspaces i we may allow for 

(Fr)^ Sr+1 = B r (a) w r e^ Sr+1 , B r (a) = e ~ 2 ^ d ^ (7) 
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where /ii.../As r +i refer to the coordinates of these subspaces and w r is an arbitrary constant. 
As for the elementary Ansatz, the sum over i in the exponent runs over the subspaces spanned 
by the form. It is easy to check that this Ansatz, already solves the form equation of motion. 

Note that in contrast to the elementary Ansatz, the solitonic field strength does not have a time 
index. Therefore the matching condition for the dimensions differs. Given, for example, a split 
(do,di) = (3,6) one has to use a 2 form instead of a three- form to fit into the 3-dimensional 
subspace. The above Ansatz then reads F^^^ = exp(— 6ao)we MlAt2A13 . 

Why did we use the terms "elementary" and "solitonic" for the above Ansatze? This is again 
motivated by a close analogy to p-brane solutions. Consider again the simple split into two sub- 
spaces, do = p + 1 and d\ = D — p — 2, where do is the analog of the p-brane worldvolume, and 
d\ is the analog of the sphere surrounding the brane. There are two types of the simple brane 
solutions. Each is a source of a form field, but the field is aligned differently in the two cases. 
For the "electrically charged" elementary brane, the form field is aligned in the hyperplane of the 
brane together with the radial direction. This has an analog in our elementary Ansatz with the 
form aligned in the do subspace together with the time direction. For the "magnetically charged" 
solitonic brane, the form field is aligned in the sphere surrounding the brane, which in turn has an 
analog in our solitonic Ansatz with the form aligned simply in the d\ subspace. 

Having specified the form of our Ansatz, we now look to solve the equations of motion derived 
from the action (Q). However, it is in fact easy to show that, under a very mild restriction on the 
structure of forms to be discussed below, the resulting equations of motion can be obtained from a 
reduced Lagrangian which depends only on «j, (f>, N and f r , each as functions of r. The Lagrangian 
is given by 



C = E 



n— 1 n—1 | 

dia? - £ did^a'j + -— — 0' 2 + V e - N 2 V S 

i=0 i,j=0 



(8) 



with 



V e = ^]>> r (a)e-^/; 2 



^B r (a)^e- fl W*+Ae- a ^ (9) 



| v-^n— 1 , 

E = —e^i=o di0li . 
N 

In the definitions of the potentials V e and V s , the sum over r is understood to run over all the 
elementary and solitonic configurations which have been chosen according to the given rules. In 
general, this even includes the possibility of having a superposition of two of these configurations 
for one form. However, in this case the nonvanishing components of the two configurations should 
differ in at least two spatial directions in order to obtain the correct equations of motion from 
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Lagrangian (||) (otherwise the equations of motion contain an additional cross term which cannot 
be derived from eq. @). We will adopt this fairly mild restriction on the structure of the forms 
thereby gaining the advantage of having a simple description of our models in terms of a reduced 
Lagrangian. 

The equations of motion for the functions f r to be derived from eq. d§[) read 

J- (EA r e- a{ - s ^ f r ) = . (10) 

The first integrals are 

f r = v r E- l A; 1 e a ^ )<t> (11) 



where v r are integration constants. Equation ( |IT| ) can be used to eliminate f r from the elementary 
potential V e . This then reduces the problem to solving the remaining equations of motion now 
given purely in terms of «j, <fi and N. 

In fact we find that the remaining equations can also be derived from a simple reduced La- 
grangian. This Lagrangian describes the motion of a particle in a "moduli space" the coordinates 
of which are the scale factors together with the dilaton. We find E is then the metric for the 
particle worldline in the moduli space. To see this structure, we first introduce the notation 
a = (aj) = {o.i,4>) for a general point in the moduli space. We also define a particular metric on 
the moduli space Gjj by 

Gij = 2(di5ij — didj) 

Gi n = G n i = (12) 
8 



G 



n it 



D-2 

The equations of motion for a and N following from eq. (|8|) then take the simple form 

±{EGcx')+E~ 1 ^- = (13) 
^Ea' T Ga' + E~ X V = 0. (14) 



U = e 2 ^=o d ^ (jLv e + V.) (15) 



The quantity U is given by 



where f r in V e has been replaced using eq. (|TT|) . Clearly, these equations of motion can be derived 
from the simple Lagrangian (it is convenient to make the change of variables from N to E) , 

C = ^Ecx' T Ga' - E~ X U (16) 



We note from its definition (12), that the metric on the moduli space has Minkowskian signature 



( — h + •••+), thus we are really free to interpret (16) as the worldline Lagrangian for particle 



moving in the moduli space. The first term is kinetic, while U defines a potential in the moduli 
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space. Further, E is the metric on the particle worldline. It encodes the gauge freedom that we 
have not yet chosen a particular time parameterization to describe the evolution of the spacetime. 



As usual its equation of motion (14) is a constraint. 

It is useful to rewrite the potential U in a more systematic way as 



U = -J^Ur exp(q r • a) (17) 
z r=l 

where the sum runs over all elementary and solitonic configurations as well as a possible cosmological 
constant term. The constants u r represent the integration constants v r in eq. (11) for elementary 
forms, the constants w r in the Ansatz (|7|) for solitonic forms or a cosmological constant. The type 
of each term is specified by the vectors q r which can be read off from eqs. (||), (|7|), (^) and the 
definition (|l5|). For an elementary 5-form they are given by 

q( el ) = (2€idi, a(5)) , ei = 0, 1 , S=J2 tidi (18) 

i=0 

with ej = 1 if the form is nonvanishing in the subspace i and ej = otherwise. For type II theories 
the dilaton couplings a(6) are given in eq. (|3|) and @). To account for the D = 11 case (or constant 
dilaton solutions) we may just set a(5) = 0. 

Let us give an example at this point. An elementary IIA RR 3 form, put into the first subspace of 
a D = 10 space split with (do, d\) = (3, 6), is characterized by a vector q = (6, 0, —1/2). It generates 
a potential term in ( |i"7| ) which depends on ao and the dilaton but not on a\. More generally, a 
potential term describing the effect of an elementary form depends only on those scale factors which 
correspond to subspaces spanned by the form. Since the entries qi of q are always positive, the 
potential tends to drive the scale factors for these subspaces to smaller values. Therefore these 
subspaces tend to be contracting or at least less generically expanding than others. 

The situation for a solitonic 5-form is in some sense complementary. It is specified by a vector 

q( so1 ) = (2etdi, -o(5)) , it = 1-^ = 0,1, ~5 = D-2-5 = J2^idi (19) 

i=0 

with ej = 1 if the form vanishes in the subspace i and ii = otherwise. For example, a solitonic 
IIB RR 2 form in the first subspace of a space split as (do, d\) = (3, 6) is specified by q = (0, 12, 1). 



It generates a potential term in (17) which depends on a,\ and the dilaton but not on ao- More 
generally, in contrast to the elementary case, the potential term now depends on those scale factors 
corresponding to subspaces not spanned by the form. 
Finally, a cosmological constant is characterized by 

q^^,-^-) . (20) 



D - 2 



Note that for all these vectors J2i=o Qi > 0, a fact which we will use later on. The moduli space 
metric allows us to define a natural scalar product on the space of vectors q.. It is the product 
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q • a that appears in the potential U. Thus q is naturally a covariant rather than a contravariant 
vector in the moduli space, so we define 

< qi,q2 >= qrG _1 q 2 (21) 



with the inverse of G given by 



(G ~ 2(D-2) + 2di Sij 



(G- l ) in = (G- l ) ni = (22) 

\S* )nn — g 

Since the metric G has Minkowskian signature, we can distinguish between space- and time-like 
vectors q. As we will see, the structure of the solutions depends crucially on this distinction. 

Generically, the models specified by the eqs. (|l3|), fll4|) and ( |i~7l ) cannot be solved. A complete 
solution, however, can be found if the potential U consists of one exponential term only or if contact 
with Toda theory can be made. These two cases will be discussed in section 3 and 4. An interesting 
observation at this stage is that the potential U may contain exponentials with different signs in 
front of the dilaton <f> and therefore may have dilaton minima at finite values. It is an important 
question how this influences the dynamics of the dilaton, in particular since the stabilization of the 
dilaton is one of the major problems of string cosmology. This will be analyzed more completely 
elsewhere [psj], 

3 Solutions with one potential term 

In this section, we will analyze models with just one form turned on (or a non-zero cosmological 
constant). The form may be elementary or fundamental and there may be any number of subspaces. 
All of these cases correspond to a potential 

U = ^u 2 exp(q • a) (23) 

where u 2 > 0. Here q is one of the vectors specified at the end of the last section. We will start by 
giving the general form of the solution and then give a simple example in section 3.2. 



One way of solving the equations of motion for a potential (23) is to use the gauge freedom in 



the definition of the time coordinate. We can always choose a gauge such that 

N = exp((d - q) • a) (24) 

where we introduce a vector giving the subspace dimensionsd = (e£j, 0). This implies E = exp(q-o;) 
and the following set of equations for a 

^-(GEa') + lu 2 q = 
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E T 1 ' ri i i 2 

-a Ga + —u 



. 



1 

2 ' 2 l 

In this form they can be integrated immediately, leading to the general solution 

ST 



where 



a = c In ti — r + w In 



2G- X q 



Ti — T 



(25) 



(26) 



(27) 



< q,q > 

The sign s = ±1 is determined by s = sign(< q, q >) Q and w, k are integration constants subject 
to the constraints 



q • w 

w T Gw 

q • k 



1 


In 



(28) 



u 



< q,q > 



T\ is a free parameter which we can take to be positive. The opposite choice t\ < leads to a 
shift in the range of r which does not affect the physics of the solution. The range of r should be 
chosen to ensure a positive argument of the second logarithm in eq. (^). This depends on the sign 
of < q, q > and we have the two cases 



(29) 



< r < n for < q, q > > 
r < or r > t\ for < q, q > < 

Which of these cases is actually realized in type II models? Using the vectors q given in the end of 
section 2, we find for a solitonic or elementary 5 form (or a cosmological constant which is similar 
to a RR (D — l)-form) 



< q, q > = 



D 



-55 



4 

D-2 



NS 
R 



(30) 



which is always positive. Also the D = 11 3 form leads to a positive result, as can be seen from the 
above formula by setting a(5) = 0. We conclude that, in the present context, we are dealing with 
spacelike vectors q only, and we have < r < t\. 

For concreteness, let us give an example of how to compute the coefficients c for a given model. 
Consider an elementary 5 form in D = 10 and a spacetime split into two subspaces with dimensions 
do = 5 and d\ = D — 5 — 1 = 5 + 1. Then from eq. (12) we find the metric 

/ -25{5-l) -26(6 + 1) \ 



G 



V 



-25{5 + l) 




-25(5 + 1) 
1 



(31) 



Here we disregarded the somewhat marginal possibility that q is a null vector, i. e. < q, q >= 0. 
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and from eq. (E^) its inverse 



/ 1 - J- 

27 16 



j_ 

16 



16 



1 



V 



2(5+1) 





16 





1 J 



(32) 



Inserting this into eq. together with the vector q = (25, 0, a(5)) characterizing the elementary 



5 form, we find the following coefficient c of the first term in the solution (26) 

_ / 8 5__ a(5) \ 

°~ {2(8 + 8)' 2(6 + 6)' 2 ) ■ 



(33) 



The above example provides useful insight into the relation of our solutions to p-brane solutions. 
As discussed in the last section, we would expect this example to correspond to an elementary (8—1)- 
brane solution. Eq. ( p6| ) has been written in a form which makes this relation transparent. The 
first term in eq. ( p6| ) is the analog of an extremal, BPS, p-brane solution. To see this, we remind 
the reader that a brane BPS solutions is characterized by a proportionality of all aj to a certain 
harmonic function h; that is a\ ~ h. The (8 — l)-brane constants of proportionality ||] are exactly 
reproduced by eq. (|33l ) . The second term in eq. (|26|) represents the analog of nonextremal extensions 
which also have been studied in the p-brane case [p4| , || |25| . The difference here is that the second 
term cannot be set to zero since wi = for all I is incompatible with the constraints (f28|). Another 
possibility is to choose t\ = which by eq. ((2^) works for timelike vectors q only. Unfortunately, 
as we saw, such vectors do not occur within D = 10 type II and D = 11 supergravity. By analogy 
with p-brane solutions one would therefore conclude that the solution (]26| ) does not preserve any of 
the supersymmetries of these theories. This can, in fact, be checked explicitly for certain examples 
by using the supersymmetry transformations of the gravitino. 

So far, our solutions have been expressed in terms of the time parameter r which is defined 
by the gauge choice (p4|) for N(t). For a discussion of the cosmological properties of our models, 
however, we should reexpress them in terms of the comoving time t, that is in the gauge where the 
N = 1. This can be done by integrating the defining relation dt = N(t)cIt. The explicit expression 
for N(t) is given by inserting the solution (^) into the gauge fixing equation (p4|) for N(r), which 
gives 

N = exp((d - q) • k)|7i - tI-^-^tI*- 1 (34) 



with 



x = d • w 

A = d.c = 2 <d ' q> 
< q,q > 



(35) 
(36) 
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Note that the quantity x depends on the specific choice of integration constants w whereas A is 
a fixed number for a given model (once the dimensions d of the subspaces and the Ansatz for the 
form encoded in q have been chosen). Depending on the values of x and A, the gauge parameter 
iV may have singularities at r = and r = t\. This determines the allowed range in the comoving 
time t as we will discuss in detail in the next subsection. 

Another quantity which is of importance in discussing the physical content of our solutions is 
the scalar curvature R. For the solution (^6|) it is given by 

fl~ |ri -r\ 2 ^- A ^T\~ 2x P 2 {w,a,T) (37) 

where P2 is a second order polynomial in r which we will not need explicitly. We have also omitted 
an unimportant constant of proportionality which depends on the integration constants. The first 
two factors in this equation indicate potential singularities at r = and r = t\, depending on x 
and A as in the case of the gauge parameter N. However, in contrast to singularities in N, such 
singularities are true coordinate independent curvature singularities. They will be further discussed 
in the next subsection. 



3.1 Cosmology of solutions with spacelike q— vectors 

As already mentioned, the case of spacelike q-vectors is the most important in our context since 
all vectors q appearing within the D = 10 type II theories and D = 11 supergravity are spacelike. 
Even under this restriction on q, our solution (|2^) covers a large number of different models. In 
this section, we will discuss the cosmological properties which can be extracted from these models 
in general. A concrete illustrating example will be given in the next subsection. 

Recall that the singularity structure of the solution (^) is determined by the quantities x and 
A defined in eq. pj3) and (B§). What values are actually allowed for these quantities? From 



< d, q >= - J2j=o Qj/ 2 i D ~ 2 ) an d E"=o Qj > [t follows from eq. © that A < 0. Though we 
are not able to give a general proof, the parameter x, which unlike A depends on the parameters 
of the solution, turns out to be either x < A or x > in all specific examples we considered. This 
divides the set of initial conditions into two disconnected subsets corresponding to two classes of 
solutions with different properties. 

We begin our discussion of these properties by analyzing the allowed ranges in comoving time 
t. Recall from eq. (p^) that the time parameter r is always in the range < r < t% since we have 



< q, q > > 0. The singularity structure of the gauge parameter N in eq. (34) then shows that this 
range is mapped to the following ranges in t 

. ,—00, iil for x < A < , (— ) branch , „ 

t£ { 1 1 1 ' . (38) 

[to, +00] for x > > A , (+) branch 



13 



Here to an d ti are t wo finite unrelated values that appear as integration constants from integrating 
dt = N(r)dr. Thus we have found two disconnected branches corresponding to asymptotically 
positive and negative time ranges. 

Let us next discuss the scalar curvature in each branch. We start with the (— ) branch. As 
inspection of eq. (|37| ) shows, the curvature vanishes as t — > — oo (r ~ 0) since x < A < 0. With 
increasing time R grows and, finally, the system runs into the curvature singularity at t = t\ 
(r = n) since the power 2{x — A) of the first term in eq. (|37j) is negative. Therefore, classically 
the (— ) solution cannot be continued beyond this point. 

In the (+) branch the situation is similar but reversed in time. At t = to (r = 0) we find a 
singularity since x > in this branch. The solution cannot be extended into the past. As t — > oo 
(t — tl) the curvature behaves smoothly and approaches zero. 

Though generically correct, the above argument has a loophole. For very specific values of the 
initial parameters w, the polynomial Pi is proportional to r or T\ — r so that it can cancel against 
one of the first two factors in eq. ([371) which cause the singularity. If \x\ is sufficiently small, the 
singularity may disappear completely. For the (— ) branch this is realized if w n = c n and x > —1/2. 
For the (+) branch it occurs if w n = and x < 1/2. This phenomenon is quite similar to what 
happens in the curvature singularity free WZW model of ref. p9| 



So far, we have considered quantities which provide information about the behaviour of the 
total D dimensional space only. Let us now turn to the individual subspaces of dimension di. To 
analyze their behaviour, we should calculate their respective Hubble parameters Hi in terms of 
the comoving time. In fact, it is possible to explicitly express the comoving time t(r) in terms 
of hyper geometric functions. It is, however, more instructive to look at the asymptotic regions 
r ~ (corresponding to t — > —00 for the (— ) branch and t ~ to f° r the (+) branch) and r ~ t\ 
(corresponding to t ~ t\ for the (— ) branch and t — > 00 for the (+) branch). In these regions the 
Hubble parameters can be written as f\ 

Hs6, = ih < 39 > 

with the constant expansion coefficients p satisfying 

pGp = , d • p = 1 . (40) 

The time shift t s depends on the asymptotic region and branch under consideration. The sign of 
t — t s , however, is always well defined : It is negative in the (— ) branch and positive in the (+) 
branch. If we combine the two equations ( f40|) and use the explicit form of the metric G in (|l2|) we 
find 

4 n ~ 1 

■pz^pl + E d ^ = 1 • ( 41 ) 

4=0 



2 The dot denotes the derivative with respect to the comoving time t. 
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The explicit expressions for p in terms of the integration constants are 

? at r ~ 



(42) 

5=A ' - T l 



They have been calculated using the general solution (|2q) and the asymptotic limits of N(r) to be 



read off from eq. (|34|), The behavior of the Hubble parameters (39) along with eq. ( f40| ) indicates 
that the solutions behave asymptotically like those of pure Kaluza-Klein theory with a dilaton. 
This can be seen by a comparison with the solutions of ref. ||. Therefore, one expects that 
the potential U provided by the form is effectively turned off in these limits. In fact, inserting 



the general solution (26) into the potential (23), we find U ~ (ri — t)t which implies that U is 
effectively zero near r ~ and r ~ t\. The effect of the form is therefore to generate a mapping 
p(r ~ 0) — > p(r ~ ri) between two pure Kaluza-Klein states. 

What do the above results mean for the evolution of the subspaces? We consider the (+) branch 
first. Remember that t — t s > in this branch so that from eq. (|3^) a positive pi results in expansion 
and a negative pi in contraction. Moreover, the equation d • p = 1 shows that at least one of the pi 
has to be positive. Consequently, at least one of the subspaces has to be expanding. From eq. ( [41~| ) 
we conclude that \pi\ < 1 always. The expansion is therefore subluminal, i. e. the scale factor e ai 
grows less quickly than the horizon size H^. This behaviour is similar to a radiation or matter 
dominated universe corresponding to pi = 1/2 and pi = 2/3, respectively. 

The situation is completely different in the (— ) branch. Since t — t s < 0, a positive pi results 
in contraction and a negative pi in expansion. Now we conclude from d • p = 1 that at least 
one subspace must be contracting. Since we are in the negative time range, eq. ( [39] ) shows that 
expansion {Hi > 0) goes along with an increasing Hi, that is, a shrinking horizon size. Such 
a behaviour is also called superinflation since scales are stretched across the horizon even more 
rapidly than in "ordinary" inflation where the horizon size is approximately constant. It is very 
similar to the superinflating expansion in pre-big-bang models |21]. However, it should be stressed 



that it was obtained in the Einstein frame as opposed to the string frame used in those models. 

Our solutions allow various patterns of expanding and contracting spatial subspaces. The details 
of the evolution depend on the partition {di}, the form and the subspaces it occupies and the initial 
conditions. Examples with 3 expanding and 6 contracting spatial dimensions as t — > oo exist, as 
were given in ref. ||. Though this is by no means preferred in the present context, its existence 
is not entirely trivial since some models fail to reproduce such a "desired" situation. The effect 
of the form can be quite dramatic. For example, it can reverse expansion and contraction of two 
subspaces during the early asymptotic period into its converse during the late period. 

A "preferred" cosmological scenario suggested by these solutions consists of a combination of 
the (— ) and the (+) branch to account for inflation as well as for a postinflationary subluminal 
expansion. It is remarkable that such a pre-big-bang scenario can be realized in the Einstein 
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frame. The apparent shortcoming of this scenario is that the (— ) and the (+) branch constitute 
two different a priori unrelated solutions. As for string frame pre-big-bang models, one might 



argue [21] that scale factor (T) duality between the branches provides the correct transition and 



some hints for this have been found in a two dimensional toy model [p(J . 
3.2 An example 

Up to this point our discussion has been rather general. Let us now illustrate the steps in our 
solution by giving a simple example. It is not meant to represent the most realistic case in the 
cosmological sense, but rather to explain our general results as concretely as possible. 

We consider the following situation : 10-dimensional spacetime is split into two subspaces with 
d = (do,di,0) = (3,6,0) and an elementary IIA RR 3 form occupies the 3-dimensional subspace. 
This implies the Ansatz 

ds 2 = -N 2 (r)dT 2 + e 2a °dx. 2 + e 2ai dx 2 



in accordance with the eqs. (|5]), (^). The equations of motion for this example can be derived from 
the Lagrangian 



C = E 

with the elementary potential V e and E given by 



6ag 2 — 30a'± 2 — 36a' a'i + -^(f)' 2 + V e 



(44) 



2 N 

which come from the general equations (||) and (||). Note that this Lagrangian depends on the 
function /' which appears in the Ansatz for the form, but not on / itself. Therefore the equation 
of motion for / can be integrated to give the first integral 

/' = u J B~ 1 e 6Q0+ ^ (46) 



with an integration constant u. From the Lagrangian ( f44[ ) we can compute the equations of motion 
for a > a i an d 4>- Using eq. Q46| ) to replace /' in these equations, we arrive at 

^(E(-12a' -36a' 1 ))+3u 2 E- 1 e 6ao - <t, / 2 = 

— (E(-36a' Q - 60a[)) = (47) 
ar 

±f E ^)--u 2 E- 1 e ea °-^ 2 = 0. 
dr v ' 4 
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Let us compare these equations with the general ones given in the moduli space formalism in (13) 



(14) and fll7|). We see that they can be indeed written in this form if we set 

/ 

G = 



\ 



-12 
-36 




-36 \ 
-60 
I) 



(48) 



and 



U 



1 



-u e 



(49) 

The matrix G above is consistent with the general formula (^) with do = 3, d\ = 6 and D = 10. 
In eq. (]l7|) we introduced a systematic way of writing the effective potential by introducing a 



characteristic vector q r for each form. From eq. (49) we read off that this vector is given by 
q = (6,0,-1/2) for our example. This coincides with what one gets by applying the general 
rule (g) to the breakup (do,d±) = (3,6) and a 5 = 3 form in the 3-dimensional subspace. The 
dilaton coupling a(S) for a RR 3-form needed in eq. (0) follows from eq. (||) to be a(5) = —1/2. 

In section 2 we also defined a scalar product (^) on the space of vectors q using the inverse of 
G. From eq. ([4S|) G^ 1 is given by 



G~ 



/ 


5 


1 





\ 




18 


16 




1 


1 









16 


18 




\ 








1 


/ 



(50) 



which agrees with the general formula ( |22|) for (do,d±) = (3,6) and D = 10. One can easily verify 
that < q, q >= qG _1 q = 4. Therefore q is indeed a spacelike vector, in agreement with the general 
result (pQ) which showed that this is true for all vectors obtained from type II forms. 



The main problem in solving the system of equations ( f47|) is the existence of two different 
exponentials, one coming from E, eq. (|45|), the other coming from the effective potential U, eq. (fi"9"|). 
Fortunately, we have a gauge freedom (time reparameterization invariance) encoded in N which we 
can use to get rid of one of the exponentials. In the next section we will use this freedom to set 
E = 1. Here, however, we choose the another possibility, namely to gauge away the potential by 
setting E = exp(6ao — 0/2). Given the definition of E in eq. (^5|), this implies 

(51) 



N = exp(-3a + 6a 2 + <t>/2) 



in accordance with the general formula (p4|) for d = (3, 6, 0) and q = (6, 0, —1/2). With this choice, 
the equations of motion (|47| ) turn into 

2 



d (e 6Ql -^ 2 (2a / 1 + 64)) = "" 



d,T 

d_ 

Ik 



(52) 
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dr V r / 2 

e 6a 1 -0/2 (6a /2 +36a , a , 2+3Oa ,2_ 0/2) = ^_ 



This is consistent with the general form ( |2q ) found for models with one term in the potential. 
Taking appropriate linear combinations of the first three equations we can derive an equation for 
the remaining exponent 6ao — </>/2, which can be solved. Then oto, ai, <j> can be expressed in terms 
of this solution. In this way one arrives at a general solution of the form ( |26| ) with coefficients c 
given by 

.=(444). ( 5 3) 



,16 16 4 

and the following constraints on the integration constants 

6wi - -w 3 = 1 
12wj + 72 Wl w 2 + 60w| = vol (54) 
6fci - -k 3 = \n(u 2 ) . 

One may arrive at the same result by inserting q = (6,0,-1/2) and G^ 1 from eq. (|50|) into the 
general formulae (pTj), (p8|). Recall that the time parameter r is restricted by < r < t\. 

To discuss the cosmology of these solutions we must perform a transformation to comoving time 
t. To do this, we need the explicit form of the gauge parameter N which we find by inserting the 



solution (gg) with fl53|), fl5j) into eq. (fID 

N = e _3fc o+ 6fc i+ fc 2/2| Tl _ r |-x+A-li ix-1 

Here x = 3wo + and A = —3/16. The quantities x, A have been generally defined in eq. (j35"l) , 
( p6] ) and their values can be easily reproduced by inserting d = (3,6,0) and c from eq. (^). The 
range of comoving time obtained by integrating dt = N{r)dr over < r < t\ crucially depends on 
the singularities in N. Eq. (55) shows that there are potential singularities at r = and r = t\. 



Their appearance is controlled by the value of x. 

Let us therefore analyze which values of x are allowed. The first two constraints (|54|) may be 
solved to express, for example, wq and w\ as a function of u>2- This shows that x depends on one 
free parameter only. Furthermore, since the second constraint ( |54"|) is quadratic in wj, we find two 
branches satisfying x < A = —3/16 and x > 0, respectively. We refer to these two branches as the 
(— ) and the (+) branch. From eq. (^) we see that < r < t\ is indeed mapped to the comoving 
time ranges given in eq. (p9|); that is to [— cx),ii] for the (— ) branch and to [to 5 oo] for the (+) 
branch (io>£i are integration constants). Moreover, the scalar curvature (|37|) has a future timelike 
singularity in the (— ) branch and a past timelike singularity in the (— ) branch. Both types of 
solutions are therefore not extendible. 



18 



Information about the evolution of the two subspaces and the dilaton can be obtained form 
the respective Hubble parameters H = a written as a function of comoving time. They can be 
calculated if r ~ or r ~ t\ since N in eq. (|5q ) can be integrated in these limits. Doing this for our 
example by using eq. (|26|), (|53|), ( pi ) and iV, r(t) calculated from eq. (plf), we find H to be of the 
Kaluza-Klein form (|39|), fl4Q). The expansion coefficients p depend on the integration constants w 
as in eq. @). 

In fact, using the first two constraints ( p4| ) we may rewrite p as a function of only, as we did 
for the parameter x before. The asymptotic expansion properties of our example therefore depend 
on one free parameter only. Instead of giving the explicit formulae, which are not particularly 
enlightening, let us give a graphical representation of p = p(w2)- We concentrate on the (+) branch 
(the expansion coefficients in the (— ) branch can be worked out analogously) where to < t < oo 
and the asymptotic regions are characterized by t ~ to ( T — 0) and t —* oo (r ~ t\). The results 
are given in fig. 1 (t ~ to) and fig. 2 (t — > oo). In both figures \po\, \pi\ < 1 always, which illustrates 
our general result that the expansion in the (+) branch is always subluminal (cf. eq. (fIT|)). 



1.5- 

1- 


Pc^_ 

V Po 


i — i — i i i — i — i n 




-2 -1 ^^0/ 

-0 75- 
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-1.5- 
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Fig 1: Expansion coefficients for the (+) branch at t ~ to- 
We see that an early expansion of the 3-dimensional subspace (po > in fig. 1) is turned 
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into a contraction as t — > oo (jjq < in fig. 2) for a wide range in W2- This can be un- 
derstood from the «o dependence of the effective potential (49). Moreover, the 6-dimensional 
space is always expanding as t — > oo (p\ > in fig. 2). In a more realistic model, such an 
expansion should be stopped by, for instance, a nonperturbative potential for the modulus a%. 




-1.5 1 



Fig 2: Expansion coefficients for the (+) branch at t —* oo. 



4 More general potentials — Toda models 

So far, we have explicitly solved the case with just one form turned on. As already mentioned, 
models with two or more forms are only soluble if they correspond to Toda models. In this section we 
will establish the relation to Toda models, and discuss two classes of solutions with some examples. 

One recalls that the main result of section 2 was that the problem of finding cosmological 
solutions could be reduced to solving for the motion of a particle in a moduli space with a particular 
potential. The coordinates in the moduli space were an n + 1 component vector a = (aj,</>) 
containing the scale factors aj, i = 0, n — 1 of the n maximally symmetric subspaces in which 
we divided the D-dimensional space and the dilaton. The dynamics of a were then described by 
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the particle equations of motion ( |13|) , fll4|) with the potential given by ( |17| ) . The inverse metric on 
the moduli space ( ]22| ) defined a natural scalar product < , > ( |2~l| ) for covariant vectors living in 
the moduli space (such as q). We also recall that the metric had Minkowskian signature. To make 
contact with Toda models we would like to write all the moduli-space equations of motion in terms 
of the scalar product. To this end we introduce a "lowered-index" covariant vector corresponding 
to the contravariant coordinate vector a by 

P = Ga. (56) 

Lets us also choose to work in the gauge 

N = exp(< d,/3 >) . (57) 

which implies that E = 1. In this gauge, the moduli space Lagrangian ( |l~6| ) written in terms 
of (3 reduces to a form Cq. The constraint ( |l4|) can then be expressed as the vanishing of the 
corresponding Hamiltonian TCq so that the whole system is described by 

Co = \<(3\(3'>-U 

Ho = j < P'iP' > + u = ( 58 ) 



U = ^E n rexp(<q r ,/3>) . (59) 

1 r=l 

The prime denotes the derivative with respect to the time parameter r. Recall that each of the 
m terms in the potential U corresponds to a certain Ansatz for a form specified by the respective 
vector q r , r = 1, m. 



A general solution for this system can be found if it corresponds to a Toda model |31[| ; that is, 
if the matrix 

<q r ,q s > 

A rs = (60) 

Q 

can be identified with the Cartan matrix of a semi-simple Lie group G for some constant q f^. Note 
that A is an m x m matrix, where m is the number of forms excited, each specified by a different 
q r . 

Let us briefly show that this property does not depend on the particular coordinate system 
in moduli space we have chosen. Under a linear coordinate transformation a — > a = Pa, the 
covariant vectors (3 and q r transform as (3 — > /3 = P~ lT (3 and q r — ► q r = P~ lT q r , while the new 
inverse metric is given by G^ 1 = PG~ 1 P T . Clearly all scalar products are unchanged. Thus, the 
matrix A rs is invariant under such an operation and the property of being (or not being) a Toda 
model is coordinate independent. In particular, the transition from the Einstein to the string frame 
3 If the group is semisimple q can be different for each simple factor. 
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discussed in the beginning of section 2 can be described by such a change of coordinates in moduli 
space. This shows that our decision to work in the Einstein frame does not affect the essentials 
of the present discussion. A similar remark applies to duality transformations within the class of 
models considered here, which we also expect to act as linear transformations on (3. 

Are there really such Toda models among the examples one obtains from type II theories? To 
answer this question systematically, and to give a simple criterion to decide in each specific case, 
we will now compute the possible entries of the matrix A. We consider two forms of degree <5i and 
5 2 with vectors qi, q2 as defined in eq. (||) and (^). From eq. (|2~TD, their scalar product is given by 



(61) 



(el.NS) fel.NS) 


= 25 12 


< q ru ci ' R) > 


D-2 

= 2S 12 Si 5 2 + 


<q (el,R) )q (el I NS) > 


= 25 u - 2 


<qi SOl U SOl) > 


= < q 5 el U el) >-2 


<qi Cl U SOl) > 


= -<qi el U el) > 



and 

(62) 
(63) 

where 5\ 2 = ^27=0 € i is the "spatial overlap" of the two forms, that is the dimension of the 

subspace in which both forms are non-zero. As before, the cosmological constant of massive IIA 
is analogous to a RR (D — l)-form. The above expressions show that the entries of the matrix A 
are small integer numbers (at least for D = 10) of either sign which are just the right properties 
to construct a Cartan matrix. For a number of models A can indeed be identified with a Cartan 
matrix. Explicit examples of this have been discussed in ref. [Q. Here we are mainly interested in 
the general structure of the solutions. 

4.1 Orthogonal vectors q 



A simple class of soluble models [32] is characterized by orthogonal vectors q r ; that is, < q r , q s >= 
for r 7^ s. Then the matrix A is diagonal and can, via appropriate rescalings, be identified with 
the Cartan matrix of SU(2) m . 

Before we solve the equations of motion, let us show that such a situation can indeed occur 
within type II theories. We consider two forms of type IIA, namely an elementary 3-form specified 
by qi and a solitonic 1-form specified by q2 occupying different subspaces. In the language of 
eq. ( |6l| ) that means S± = 3, 5 2 = 1 and 5± 2 = (no overlap). Then from eq. (|63|) and ([H]) we have 
< qi,qi >=< q2,q2 >= 4 and < qi,q2 >= and, hence, we have orthogonal vectors q r . 

Let us now solve the equations of motion. It is useful to introduce an orthonormal basis e,, 
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i = 0, ...,n in moduli space satisfying 



< ei,ej >= r)i5ij (64) 

with r/Q = — 1 and iji = 1 for i > 0. In the following, we will use indices ... = 0, ...,n to label 
these basis vectors. Remember that the scalar product has a Minkowskian signature so that one 
of these vectors, eo, is normalized to —1. We have assumed that the vectors q r , r = 1, ...,m are 
orthogonal. This allows us to identify some of the basis vectors e^, with a normalized version of 
the q r 

e r = — — , q r = \J\< Qr,qr > I , r = l,...,m. (65) 

Qr 

Note here, that the q r are orthogonal spacelike vectors. This implies that their number (the number 
of forms) is less than the dimension of the moduli space; that is, m < n. Furthermore, none of 
the q r can be identified with the timelike basis vector eo- The basis vectors ( |64l) therefore fall into 
two groups, those which can be identified with forms and those which cannot. Correspondingly, 
we partition the index % = 0, n as i = (a, r), where we use the same index r = 1, m as for the 
forms to label the basis vectors proportional to q r and the index a = 0, m + 1, ...,n to label the 
other basis vectors. 

We expand the "dynamical" vector (3 in terms of the basis 

n 

P = Y,Pi e i- ( 66 ) 

with (time dependent) expansion coefficients Pi(r). Inserting this expansion into C$ given in 
eq. ( |5"8| ) , we find a complete decoupling in the eigenmodes pi and the resulting equations of motion 

Pr = \qru 2 r e~^ (67) 

Pa = 0. (68) 

Their solution given by 

Pr = q- l \n(g r ) 

Pa = k a {r-T a ) 



with where we define 



2 . // .. 1 i 



g r = K r cosh (z r ) , K r = , z r = -\k r \q r (r - r r ) , (70) 

Here r runs over the full real axis and k{ and Tj are integration constants. Bringing everything 
together, the solution for the original vector a reads 



a = J2piG- 1 e i (71) 

i=0 



with pi as indicated above. 



23 



The application to an explicit example is straightforward. First one should calculate the matrix 
G _1 from eq. (p2j). Then one determines the vectors q r specifying the forms from eq. (||), (|7|) 
and calculates the basis vectors e r , r = l,...,m via eq. (|65|). Finally, one extends this subset of 
basis vectors to a full orthogonal basis of the moduli space satisfying eq. (|4|) . All these quantities 
inserted into eq. ( |7ll) provide the explicit solution. 



So far we have not considered the Hamiltonian constraint TLq = in eq. (|5q). For the above 
solution it leads to 

1 n 

Ho = ? E^ ? = - ( 72 ) 

£ i=0 

Correspondingly, the constants k{ can be interpreted as the "energy" contribution of the modes p^. 
To understand the meaning of the "time shifts" r r it is instructive to compute the potential U for 
the above solution 

-. m -, m 

r=l r=l 

It indicates that the form r is operative around r ~ r r and switched off otherwise so that we 
expect asymptotic Kaluza-Klein behavior at r — > ±oo. Suppose we time order the time shifts 
T r ; that is, r n < t T2 < ... < r Tm . Then the system starting out at r — > — oo will go through a 
sequence of "transformations" each induced by one of the forms before it reaches the asymptotic 
region at r — » oo. If the time differences T ri+1 — T Ti are all of order one, then the forms will act 
almost simultaneously. However, if there is a long interval between two of the time shifts, for 
example r r2 — r n S> 1, the potential is effectively turned off during much of the intermediate period 
r n < t < t r2 . Depending on the number of long intervals, we may have up to m — 1 of these 
intermediate Kaluza-Klein regions. 

Though our solutions cannot be generally expressed in terms of the comoving time in closed 
form, this can be done in the Kaluza-Klein regions. For the gauge parameter N we find from 
eq. (j57|) and (§§ 

N = exp A lP ^j , A; =< d, ej > . (74) 
Using the explicit expression for pi, eq. (39), this can be approximated in any one of the Kaluza- 



Klein regions as 

N~N <r-, i/ = ^A; a A a + ^(±)|A ;r |A r . (75) 

a r 

Note that there is a sign freedom in the second term of the second equation. This freedom serves 
to specify the asymptotic region for which N is computed. Let us concentrate on the asymptotic 
region specified by r Ti < r < T ri+1 . Then all terms with r = ri,...,rj get a + sign and all terms 
with r = rj+i, r m get a — sign. More intuitively, a term r gets a + sign if the rth form 
has been "passed" and a — sign otherwise. Of particular importance are the two parameters v± 
corresponding to r — > ±oo which determine the relevant range of the comoving time t. We get v + 
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from the above formula for u by taking all signs to be + and u- by taking all signs to be — . From 
eq. (frB]), the quantities v± do not seem to be restricted in range. We have, however, to consider 
the Hamiltonian constraint ( |72|) on the integration constants k{. Using this constraint, we find for 
all explicit examples that u + and v_ always have the same sign. Then integration of dt = N{r)dr 



with N as in eq. (75) leads to two unconnected branches with 

u± < (— ) branch 




(76) 

u± > (+) branch 



The Hubble parameters in each Kaluza-Klein region take the standard form (39), with the expan- 
sion coefficients p expressed in terms of the integration constants as 



(77) 



The signs in the second term have to be chosen according to the same rule as above. Depending 
on the choice of these signs, eq. ( |77j ) specifies the expansion parameters in each of the up to m + 1 
asymptotic regions. To get this result we have used the definition of H in eq. (|39|) , the solution (fri|) 
and N, i(r) to be calculated from eq. (|75|) . Clearly, these coefficients fulfill the two equations (|40|) 
characteristic for a Kaluza-Klein region. Therefore, the discussion of cosmological properties given 
in section 3 based on these equations applies here as well. In particular, the expansion in the (— ) 
branch is superinflating for Pi < and the expansion in the (+) branch is always subluminal since 
\Pi\ < 1- 

In summary, we have found the direct generalization of the simple case with one form : the 
solutions split into two different branches (+), (— ) with properties similar to those discussed in 
section 3. We also have two Kaluza-Klein regions, one at early and one at late time, in each 
respective branch. If all time shifts r r are of the same order, these two regions are the only Kaluza- 
Klein regions, as in the case of one form. The transition, however, is then generated by the combined 
effect of all forms. If there are long intervals between the r r we may have up to m — 1 intermediate 
Kaluza-Klein regions, where m is the number of forms. In each region the forms are effectively 
turned off. A transition is generated when one or more of the forms become operative. 



4.2 A Toda example 

Let us illustrate the general results above by an explicit example. We break up the space as 
d = (3, 2, 4, 0) and consider an elementary IIA 3 form in the do = 3 subspace and a solitonic 
IIA 1 form in the d\ = 2 subspace. From the eqs. (^) and (^) we find the corresponding vectors 
qi = (6, 0, 0, —1/2) for the 3 form and q2 = (6, 0, 8, 3/2) for the 1 form. In the notation of eq. ( |Bl"D 
we therefore have <5i = 3, 82 = 1 and 612 = 0. Inserting these quantities into the eqs. (|6l|) , (|62|) and 
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we find < qi,qi >=< q2,q2 >= 4 and < qi,q2 >= 0. This shows that our example indeed 
corresponds to an orthogonal model. By using the inverse matrix 



G 



-i 



V 



5 


1 


1 


48 


16 


16 




3 


1 


_ T6 


16 


16 


1 


1 


1 


16 


16 


16 












\ 



1 / 



(78) 



obtained from eq. (22) this can also be verified explicitly. Apart from G~ l we need the set ej, 
i = 0, 1,2,3 of orthogonal vectors in moduli space to get the explicit solution from eq. (|7l|). The 
two basis vectors ei, e2 in the direction of the forms follow easily from eq. ( |6"5| ) by normalizing 
qi and q2. The two remaining vectors eo, e3 should be chosen to complete this to a basis of the 
moduli space. A possible choice for the complete basis is 



e 
ei 

e?, 



(6,2,8,0) 
-3,0,0,- 



-3,0, 



(79) 



A/6(l,0,2,-i 



It is straightforward to verify that this set of vectors indeed fulfills the normalization condition ( |64| ) . 
Inserting G~ l and e« into eq. (|7T[), along with the functions pi, gives the explicit solution for our 
example. It follows from eq. ( |69|) that the functions pi, i = 0, 1, 2, 3 read 

ln(g r ) 



for r = 1, 2 where g r is given in eq. (|70|), and 

Pa = k a (r- T a ) 

for a = 0, 3. The effective potential U, eq. (|73|), specializes to 

U = ~(kf cosh~ 2 (zi) + k\ cosh _2 (z 2 ) 



(80) 



(81) 



(82) 



for our example, where z r has been defined in eq. (70). Note that the first term in this potential is 
related to the 3 form and the second one to the 1 form. 

For a discussion of cosmological properties the numbers A, defined in eq. ([74]) are needed (see 
for example eq. (|75[)). From the above basis vectors they are easily computed to be Ao = — 1, 
Ai = 3/16, A2 = 7/16 and A3 = — 3\/6/16. Now we have explicitly given all quantities required 
to compute the gauge parameter N, eq. fl75|), and the expansion parameters p, eq. ([77]), as a 
function of the integration constants fcj. Note that the Hamiltonian constraint (]72l ) should be used 
to eliminate one of these integration constants. 
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Let us discuss the various cases that can occur within our model. Recall first that each form 
r = 1,2 is operative at r ~ r r only, where r r are free integration constants. The dashed line in 
fig. 3 shows the potential ( p^ ) for t\ = —1, t 2 = 1, so that the difference in the time shifts is of 
order one. It can be seen that the effect of the two forms cannot be separated in time in this case. 
We have only two asymptotic regions, r — > — oo and r — ► oo. As an example, we want to compute 
the quantities v±, which specify the behaviour of the gauge parameter N, in these regions. Using 
the constants Aj given above we find from eq. (f75[) 



3\^6 3 7 

v- = -ko —h - —\h\ - —\k 2 \ (83) 

lb lb lb 

for the asymptotic region r —* — oo. Correspondingly, for the region r — > oo we find 

, 3V6. 3 . . 7 , , nA , 

v+ = -fa - —— k 3 + — \ki + — \k 2 . 84 
lb lb lb 

Note that the r = 1, 2 terms are both negative for n_ and both positive for n + . This is in agreement 

with the general description of the sign choice, which we have given below eq. ([75]). The Hamiltonian 
constraint ([7^) reads 

H = -kl + kl + kl + kl = Q (85) 

and can be used to eliminate, for example, ko from the above expressions for v±. Doing this, it 
can be seen explicitly, that v+ and f_ always have the same sign, independent on the values of 
ki,k2,ks. As stated before, this leads to the two branches indicated in eq. (l76|). The computation 
of the expansion coefficients p in eq. (|77| ) can be carried out in a similar way. We will, however, 
not give the explicit formulae here. 

If \t2 — T\\ S> 1 another intermediate Kaluza-Klein region appears. This is illustrated by the 



solid line in fig. 3 where we have plotted the potential U, eq. (73), for t\ = —10 and t 2 = 10. 
For this choice of r r the system in the intermediate region has already been influenced by the 3 
form but not yet by the 1 form (of course one can have the converse situation by taking t\ S> t 2 ). 
Correspondingly, for the expansion coefficients in eq. ( |77| ) one should take the + sign for the r = 1 
term and the — sign for the r = 2 term. The expression for v in this intermediate region therefore 
reads 

3\/6 3 . . 7 . . . . 

^ = -ko ~ -r^k 3 + —\ki\ - —\k 2 \ . (8b) 
lb lb lb 

So far, we have carried out the discussion in terms of the time parameter r. Upon mapping 
to the comoving time t via eq. (^fj), we have two different branches, the (+) branch and the (— ) 
branch. The above discussion therefore applies to both branches. In particular, in both branches 
we may have either two or three Kaluza-Klein regions, depending on the values of the time shifts t\ 
and t 2 . Though we have not computed the expansion coefficients p explicitly, we know from their 
general properties MQ) that the evolution in each of these Kaluza-Klein regions will be similar to 
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what we have described in section 3 for one form only. 
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Fig 3: Effective potential U as a function of the time parameter r for the example described in the 
text. The solid line corresponds to t\ = —10, t% = 10, the dashed line to t\ = —1, t<i = 1. For 
simplicity we have taken k\ = hi = 1. 



4.3 Nonorthogonal vectors q 

We now discuss the general case where the vectors are not orthogonal. As in the previous 
subsection we start with an expansion (^) for (3 in terms of a basis ej, i = 0, n in moduli space. 
Again, we would like to identify some of the basis vectors with the vectors q r describing the forms. 
Therefore we define a basis with the following properties 

< e a ,e 6 > = rj a 5 ab 

<e a ,e r > = (87) 

e — 

q 

Note that the basis vectors e r are no longer orthogonal to each other. We remind the reader that 
the different types of indices were defined as = 0,...,n for the whole set of basis vectors, 

r, s, ... = l,...,m for the basis vectors corresponding to forms and a, b, ... = 0, m + l,...,n for the 
other basis vectors. With the above properties of ej, insertion of the expansion (36) into the 
Lagrangian ( |58| ) leads to 

£ = + ^Y, A rsp'rP's-U 

a " r,s 



u 



^^ufexph^AsPs) ■ (89) 
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Recall, that A rs defined in eq. (60) is assumed to be the Cartan matrix of a semi-simple Lie 
group G. Prom the above Lagrangian, we see that by introducing the basis e, we have achieved a 
separation between the trivial modes which do not occur in the potential U and the nontrivial ones 
in the directions of the vectors q r . 

The solution of the equations of motion to be derived from Cq can be written in the form 

Pa = k a (T-T a ) 

p r = -2A rr ~ l ln(g r ) 

with g r given by 

g r = b rW ex P( A • kr — A • t) . (91) 
AeAr 

Here A r are the weights of the fundamental representations of the group G, k = (k r ) is a constant 
vector restricted to the Weyl chamber of G and the set of time shifts r = (r r ) is arbitrary As before 
the constants k r are subject to the Hamiltonian constraint. The coefficients b r (X) are determined 



in terms of k r and properties of the group G but the general expression [31] is not very helpful. In 
practice, one may just insert the above solution into 



Pr + ex P A rsPs^ = 



(92) 



to find b r (X). 



Let us proceed with the simplest nontrivial example G = SU(3), that is, with a model having 
two forms excited and where the matrix A takes the form f3£ 




A = ^ (93) 



How can such a model be realized within D = 10 type II theories? We choose an elementary IIA 
3-form with vector qi and an elementary IIA NS 2-form with vector q2 such that they do not 



overlap in the spatial directions. In the language of eq. (SI) we then have 8\ = 3, 62 = 2 and 
612 = which results in < qi,qi > = < q2,q2 >= 4 and < qi,q2 >= —2. The relation ( |60| ) is 
therefore indeed fulfilled for the above SU(3) Cartan matrix if we choose q = 2. 

The weight systems of the two fundamental SU (3) representations are Ai = {(1, 0), (—1, 1), (0, — 1)} 
for 3 and A2 = {(0, 1), (1, — 1), (—1,0)} for 3. The constants k = (^1,^2) are constrained by 
2k± — k,2 > 0, 2/c2 — k± > (k is in the Weyl chamber) which ensures that the argument g r of the 
logarithm in eq. (|90[) is positive. From eq. (|92]) we get for the coefficients b r (X) 



hO-,0) = qul^p^ 62(0,1) = qu\^^ 

= qu\^ 62(1,-1) = qu\^ (94) 
6l(0,-l) = qul 2 -^ 6 2 (-l,0) = qul 2 -^ 
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with P = (2/c2 — k\){k\ — k2){2k\ — A^). The Hamiltonian constraint 71$ = in eq. (58) turns into 

1 

2 



Ko = lY, r lakZ + l(kl-k 1 k2 + k%)=0. (95) 



a 

by inserting the complete solution. 

Though significantly more complicated, the structure of these solutions is similar to what we 
found in the case of orthogonal vectors q. For r —* ±00 we have two asymptotic Kaluza-Klein 
regions. Having chosen appropriate constants k and r one may also have intermediate Kaluza- 
Klein regions. The transition between these regions, however, is more complicated and the effect 
of the two forms can no longer be separated in time. 

5 Including curvature 

Thus far, we have assumed that the maximally symmetric subspaces into which we have split the 
total space are flat. In this section, we would like to show that the more general case of curved 
subspaces can easily be incorporated into our general framework. This is of particular interest 
because it allows us to investigate a class of cosmologies, directly related to black p-brane solutions, 



some examples of which have recently been given by Behrndt and Forste 16] and Poppe and 
Schwager pC]] . 

Suppose, that instead of (B) we start with a metric 

n-l 

ds 2 = -N 2 (r)dT 2 + e 2ai dn 2 Ki , (96) 

8=0 

where d£l 2 K . is the metric of a di dimensional space with constant curvature K{ = —1,0 or +1. The 
effect of these terms is to create an additional potential V c for the scale factors on. The effective 
Lagrangian (|8|) is modified to 

C c = C - N 2 EV C , (97) 

where V c is explicitly given by 

V c = -J2 2K ie - 2a > . (98) 

i 

This potential is formally very similar to the one provided by the forms (see the eqs. @). It leads to 
terms of the form u 2 exp(q-o;) in the effective potential U which, in fact, we have already considered 
by using the general form ( |l7|) of U. All we have to do is to specify the vector characterizing a 
curvature in the kth subspace. From eq. (|98|) it follows that 



q* = (2do,...,2(d fc -l), ...,2^-1,0) . (99) 

Moreover, the constant u 2 in front of such a curvature term in U is given by u? = — 2JQ. so that it 
can be of either sign depending on the sign of the curvature K^. This is in contrast to terms in U 
describing forms which are always positive. 
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All methods presented in this paper for a potential U describing forms only can now be applied 
to the extended version including the curvature terms. The relevant information is encoded in the 
vectors q which we have specified for forms in eq. @, (0) and for curvature terms in eq. (^) above. 

For example, as in the case of flat subspaces, a model is soluble if it corresponds to a Toda 
model. To check this, one has to compute the scalar products between all the vectors q now also 
including those describing curvature in analogy to eq. (|6l|). This is easily done given the above 
expression fl99|). For instance, the scalar product between a vector q^ describing curvature in the 
subspace k and a vector q( el ) for an elementary form is given by 

<q fe ,q (cl) >=-2e k , (100) 

where we have used eq. (|2l|), ( |22| ) and eq. Recall, that = 1 if the subspace k is occupied 
by the form and e& = otherwise. Therefore, a model with one elementary form and one curved 
transverse subspace has orthogonal q vectors and can be solved following the lines of section 4.1. 
More generally, if a Toda model with curvature has been found by calculating all scalar products, 
it can be solved following the general procedure which we have described in section 4. 

By way of a specific example with a non-flat subspace, let us consider the cosmology with non- 



trivial NS three-form given by Behrndt and Forste [16]. Note that, unlike all previous examples, 
this is in five not ten dimensions. The solution is derived directly from the five-dimensional black 
hole solution with a NS three-form ]^4j by going to a regime where the time and radial coordinates 
exchange roles. The metric is given by 

dt 2 -1 + t 2 It 2 

ds2 = - (-1 + -£/*») + T^k^ + tW (101) 

where dQ 2 is the metric on a three-sphere, while the dilaton and the NS form are given by 

<t> = -- In (l - t 2 _/t 2 ) + 0o = 2M-*~%im« j (102) 

where Hi,fi2 an d /^3 are coordinates on the three-sphere. We see immediately that this solution is 



in the form of our general Ansatz. Comparing with (96), we see that the space has been split into 
two subspaces: a one-dimensional subspace parameterized by y and a three-dimensional spherical 
subspace which together with the time direction constitutes our observed spacetime. Thus we would 
write do = 1 and d\ = 3 while Kq = and K\ = 1. The NS form field spans only the three-sphere, 
without a component in the time direction. This implies a solitonic Ansatz for H^ lfl2 ^ 3 . Using eqs. 
Qi~8| ) and (|99|), we can immediately write down the vectors q# and <\k corresponding to the NS 
field and the non-zero subspace curvature, 

q H = (2,0,8/7) q x = (2,4,0). (103) 

If uh is the charge of the NS field, then, by the discussion above, the potential U for this Ansatz 
is given by 

U = ^u 2 H e tiH - a - 2e qA "- Q . (104) 
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Similarly, we can calculate the metric on moduli space and its inverse for this five-dimensional 
Ansatz using eqs. (fi^) and (p^). For the inverse metric we find 



-loo 



G~ = 

I 1/ 

This allows us to compute the scalar product between q^ and qx- We find that 



(105) 



<q//,qx>=0, (106) 

so that, in fact, we have a model with orthogonal q vectors. As a result, it corresponds to a 
simple SU(2) Toda model and we can solve the equations of motion exactly following the discussion 
given in section 4.1. Choosing the correct time parameterization, we thus re-derive the solution 
of Behrndt and Forste given above. We note that there are many other models with a spherical 
subspace which could be solved in this way. It would be interesting to know if these can all be 
directly related to black p-brane solutions. 

It is also interesting to note that the presence of a curved subspace significantly alters the form 
of the solution. We note that, unlike the previous cases, with a spherical subspace the corresponding 
term in the potential U is negative. The effect, in a solution such as that of Behrndt and Forste, is 
that there is no longer necessarily either an initial or a final singularity. From the form of the metric 
in eq. ( |101| ), it would appear that the solution is divergent at t = t + and t = t_. However, in fact, 
the scalar curvature is finite at each point (essentially since they correspond to the two horizons of 
the origin black p-brane solution). This raises the very interesting possibility that solutions of this 
form may resolve the graceful exit problem inherent in the cosmological solutions given thus far, 
providing an inflating solution which does not end in a curvature singularity. 



6 Conclusion 

Throughout the paper we have noted the close relation between the cosmological solutions we have 
been discussing and supergravity p-brane solutions. The relationship consists of exchanging the 
time coordinate for a transverse spatial coordinate, thus transforming a cosmological space where 
everything evolves in time into a brane world-volume where the fields all depend on the transverse 
radial coordinate, and vice versa. This relationship is very suggestive. In particular, it raises the 
possibility that the statistical mechanics of Hawking radiation in cosmological spacetimes (which 
is responsible for large scale structure) could have a similar string theoretical origin as the entropy 
of black holes [|J (which can be computed from the associated D-branes). We are studying this 
possibility and will report on it elsewhere. 
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Note: Shortly after this paper was finished, an interesting paper by Larsen and Wilczek [34] 
appeared which emphasizes that the solutions of Behrndt et. al. [16], Poppe et. al. [20] and some 
of those discussed in this paper actually exactly correspond to the interior of certain p-brane black 
holes, and gave several new cosmological solutions in this context. They independently point out 
that these solutions might evade the problem of curvature singularities. 
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